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Abstract 



Results are reported from a search for new physics processes in events containing a 
single isolated high-transverse-momentum lepton (electron or muon), energetic jets, 
and large missing transverse momentum. The analysis is based on a 4.98 fb -1 sam- 
ple of proton-proton collisions at a center-of-mass energy of 7 TeV, obtained with the 
CMS detector at the LHC. Three separate background estimation methods, each rely- 
ing primarily on control samples in the data, are applied to a range of signal regions, 
providing complementary approaches for estimating the background yields. The ob- 
served yields are consistent with the predicted standard model backgrounds. The 
results are interpreted in terms of limits on the parameter space for the constrained 
minimal supersymmetric extension of the standard model, as well as on cross sec- 
tions for simplified models, which provide a generic description of the production 
and decay of new particles in specific, topology based final states. 
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1 Introduction 

This paper reports results from an updated and improved search for new physics processes in 
proton-proton collisions at a center-of-mass energy of 7 TeV, focusing on the signature with a 
single isolated lepton (electron or muon), multiple energetic jets, and large missing momentum 
transverse to the beam direction The data sample was collected by the Compact Muon 
Solenoid (CMS) experiment during 2011 at the Large Hadron Collider (LHC) and corresponds 
to an integrated luminosity of 4.98 fb -1 , roughly one hundred times larger than the sample 
used for our previous search JTJ . 

The lepton + jets + ftp signature is prominent in models based on supersymmetry (SUSY) U2J- 
|7| . In R-parity-conserving models (HI, SUSY particles are produced in pairs, and their decay 
chains end with the lightest supersymmetric particle (LSP). In some scenarios, the LSP is a 
neutralino (x°), a heavy, electrically neutral, weakly interacting particle with the properties of 
a dark-matter candidate 0. The presence of two such LSPs in each SUSY event typically leads 
to a large missing transverse momentum, depending on the details of the SUSY mass spectrum. 
The isolated lepton indicates a weak decay of a heavy particle, such as a W boson or a chargino 
(X )• Multiple jets can be produced in complex decay chains of SUSY particles. This signature 
arises in many SUSY models, including the constrained minimal supersymmetric extension of 
the standard model (CMSSM) lHUl [TT| , and in simplified models |12T[l"5l , which are based on 
simplified mass spectra and decays of new particles. Both of these frameworks are used to 
interpret the results. Searches in this or similar channels have been reported by CMS |[TJ [16l 
and ATLAS llZnlll- 

Searches for SUSY particles are complicated by the presence of standard model (SM) back- 
grounds that can share many of the features of signal events. In the single-lepton final state, 
backgrounds arise primarily from the production of tt and W+jets events, with smaller con- 
tributions from Z+ jets, single-top quark production, and QCD multijet events. In the event 
topology studied here, a large observed value of f*p in a standard model event is usually gen- 
uine, resulting from the production of one or more high-momentum neutrinos. A smaller con- 
tribution to events in the high-^j tail in this search can arise from the mismeasurement of jets 
in high cross section processes such as QCD multijet events. To determine the contributions 
from these backgrounds, we use methods that are primarily based on control samples in data, 
sometimes in conjunction with specific information from simulated event samples or from ad- 
ditional measurements that provide constraints on the background processes. 

Three complementary methods are used to analyze the data, providing valuable cross-checks 
and probing different signal regions. The Lepton Spectrum (LS) method was used in the CMS 
single-lepton [1J and opposite-sign dilepton |20"| SUSY searches performed using the 2010 data 
sample. It uses the observed lepton transverse momentum (pj) spectrum and other control 
samples to predict the £t distribution associated with the dominant SM backgrounds. This 
method is sensitive to SUSY models in which the £t distribution is decoupled from the lepton 
pj spectrum, as is the case when two undetected LSPs produce a large missing transverse mo- 
mentum. The Lepton-Projection Variable (Lp) method uses the hp variable, which was developed 
for the CMS measurement of the W polarization in W+jets events [21]. This variable, described 
in Section |6j is correlated with the helicity angle of the lepton in the W-boson rest frame. Both 
the Lp and the LS methods take advantage of well-understood properties of the W polarization 
in tt and W+jets events for the background determination. The ANN method uses an artificial 
neural network discriminant built from several kinematic quantities. The ANN discriminant is 
then used in conjunction with Ifcp to define signal and sideband regions, from which the back- 
ground yield is determined. A key variable in the ANN is Mp, an approximate invariant mass 



2 



3 Data and simulated event samples 



of the system comprising the lepton and the computed with the momentum components 
transverse to the beam direction. Background events usually have Mj < M(W), where M(W) 
is the W boson mass, because the observed £x is associated with the neutrino from W — > £v 
decay. 

This paper is organized as follows. Sections [2] and [3] describe the CMS detector and the event 
samples. The event preselection requirements that are common to all methods are discussed in 
Section [4] Sections |5j |6j and [7] describe the LS, Lp, and ANN methods, respectively, for obtain- 
ing SM background estimates from control samples in data. The observed yields in data are 
compared with the background estimate obtained for each method. Systematic uncertainties 
are described in Section [HI Finally, the results, interpretation, and conclusions of the analysis 
are presented in Sections |9| and [10| 

2 The CMS detector 

The CMS detector, described in detail in Ref. |22|, is a multipurpose apparatus designed to 
study high-px physics processes in proton-proton collisions, as well as a broad range of phe- 
nomena in heavy-ion collisions. The central element of CMS is a 3.8 T superconducting solenoid, 
13 m in length and 6 m in diameter. Within the magnet are (in order of increasing distance 
from the beam pipe) high-precision silicon pixel and silicon strip detectors for charged particle 
tracking; a lead-tungstate crystal electromagnetic calorimeter for measurements of photons, 
electrons, and the electromagnetic component of jets; and a hadron calorimeter, constructed 
from scintillating tiles and brass absorbers, for jet energy measurements. Beyond the magnet is 
the muon system, comprising drift tube, cathode strip, and resistive-plate detectors interleaved 
with steel absorbers. Most of the detector systems are divided into subsystems that cover the 
central (barrel) and forward (endcap) regions. The first level of the CMS trigger consists of cus- 
tom hardware processors that use information from the calorimeter and the muon system to 
select up to 100 kHz of the most interesting events. These events are then analyzed in the High 
Level Trigger (HLT) processor farm, which uses information from all CMS detector systems to 
reduce the event rate to about 300 Hz. 

In describing the angular distribution of particles and the acceptance of the detector, we fre- 
quently make use of the pseudorapidity, rj = — ln[tan(0/2)], where the polar angle 6 of the 
particle's momentum vector is measured with respect to the z axis of the CMS coordinate sys- 
tem. The z axis points along the direction of the counterclockwise-moving proton beam; the 
azimuthal angle (p is measured in a plane perpendicular to this axis. The separation between 
two momentum vectors in tj-ip space is characterized by the quantity AR = (At]) 2 + (A(p) 2 , 
which is approximately invariant under Lorentz boosts along the z axis. 

3 Data and simulated event samples 

The data samples used in the analysis were selected using triggers based on lepton pj, and 
the transverse momenta (p T ) of the observed jets /. The overall level of jet activity was mea- 
sured with the quantity H T rigser = Ya P x , the scalar sum of jet transverse momenta satisfying 

p ! T > 40GeV. The missing transverse momentum ^ T tn SS er wa s computed in the trigger us- 
ing particle-flow algorithms 11231 124]| . To maintain an acceptable trigger rate, the thresholds on 

£ T trigger^ kpton p ^ an( j H «8S« 

were raised as the LHC luminosity increased over the course of 
the data collection period. The highest thresholds applied in the muon trigger selection were 
£ T trigger > 50 GeV ^ muon pj > 15 GeV/ and H ^ss^ > 300 GeV For e i ectron triggers, the high- 
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est thresholds applied were £ T tri ss er > 50 GeV, electron p T > 15GeV and Hj 1 gger > 250 GeV; 
a loose electron isolation requirement was also applied to help control the rate. The offline 
analysis requirements for both muon and electron events are more restrictive than those used 
in the trigger. 

The analysis procedures are designed using simulated event samples. Except for certain scans 
of the SUSY parameter space discussed later, the detector simulation is performed using the 
Geant4 package H2511 . A variety of Monte Carlo (MC) event generators are used to model the 
backgrounds. The QCD multijet samples are generated with the pythia 6.4.22 [26| MC gen- 
erator with tune Z2 [27]. The dominant background, tt, is studied with a sample generated 
using MadGraph 5.1.1.0 [28]. The W+jets and Z+jets processes are also simulated with Mad- 
GRAPH. Single-top (s-channel, i-channel, and tW) production is simulated with POWHEG [29]. 
To model the effect of multiple pp interactions per beam crossing (pileup), simulated events 
are generated with a nominal distribution of multiple vertices, then reweighted to match the 
distribution of the number of collision vertices per bunch crossing as measured in data. 

Event samples for SUSY benchmark models are generated with PYTHIA. As example CMSSM 
scenarios, we use LM3 and LM6, which are among the standard benchmarks [30] used in CMS. 
The CMSSM benchmarks are described by the universal scalar mass parameter m.Q, the univer- 
sal gaugino mass parameter mi/ 2 , the universal trilinear soft-SUSY-breaking parameter Aq, the 
ratio of the two Higgs-doublet vacuum expectation values tan /3, and the sign of the Higgs 
mixing parameter \i. The LM3 (LM6) benchmark is described by tuq = 330 GeV (85 GeV), 
m 1/2 = 240 GeV (400 GeV), A = OGeV (OGeV), tan£ = 20 (10), and y, > (0). For LM3, the 
masses of the gluino and squarks are very similar (^600 GeV), except for m (t) ~ 440 GeV, while 
the mass of the LSP is m(x%) = 94 GeV. The LM6 spectrum is heavier, with m(g) ~ 930 GeV, 
m(q) « 800 GeV, m(t) « 650 GeV, and m(jifi) « 160 GeV. The next-to-leading-order (NLO) 
cross sections for these models are approximately 4.8 pb (LM3), and 0.4 pb (LM6). 

The ANN method uses the LM0 model [30] to train the neural network. Because of its large 
cross section (54.9 pb at NLO), LM0 has already been excluded HI, but its kinematic distribu- 
tions still provide a reasonably generic description of SUSY behavior with respect to the vari- 
ables used in the neural network. The parameters for LM0 are mo = 200 GeV, mi/2 = 160 GeV, 
A Q = -400 GeV, tan $ = 10, and ji > 0. 

The results are interpreted in two ways: (i) as constraints on CMSSM parameter space and (ii) 
as constraints on cross sections for event topologies described in the framework of simplified 
models. In both cases, a large number of simulated event samples are required to scan over the 
relevant space of model parameters. For this reason, the scans are performed with the CMS 
fast simulation package [31], which reduces the time associated with the detector simulation. 

Both the LS and Lp background determination methods rely on knowledge of the W-boson 
polarization in W+jets and in tt events. The polarization effects are well modeled in simu- 
lated event samples, which are used in conjunction with control samples in data. The angular 
distribution of the (positively) charged lepton in the W + rest frame can be written as: 



where f + \, and /o denote the polarization fractions associated with the W-boson helicities 
+1, —1, and 0, respectively. The angle 0| is the polar angle of the charged lepton in the W + 
rest frame, measured with respect to a z axis that is aligned with the momentum direction of 
the W + in the top-quark rest frame. The polarization fractions thus determine the angular 
distribution of the lepton in the W rest frame and, together with the Lorentz boosts, control the 



dN 



- = f +1 | (1 + cos 9} f + /_! | (1 - cos ^) 2 +/ | sin 2 B* t , 



(1) 



dcos0 
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4 Event preselection 



px distributions of the lepton and the neutrino in the laboratory frame. 

The W polarization fractions in top-quark decays have been calculated |32| with QCD cor- 
rections to next-to-next-to-leading order (NNLO), and the polarization is predominantly lon- 
gitudinal. For t-> bW+ these fractions are f = 0.687 ± 0.005, /_i = 0.311 ± 0.005, and 
f+i = 0.0017 ± 0.0001. These precise calculations reduce the uncertainties associated with 
the W polarization in tt events to a low level. The theoretical values are consistent with mea- 
surements from ATLAS El, which obtained f Q = 0.67 ± 0.03 ± 0.06, /_ x = 0.32 ± 0.02 ± 0.03, 
and f + i = 0.01 ± 0.01 ± 0.04, expressed for the W + polarizations. 

The W polarization in W+jets events exhibits a more complex behavior than that in tf produc- 
tion. Both CMS [21 j and ATLAS [34 1 have reported measurements of these effects, which are 
consistent with ALPGEN (351 and MADGRAPH (28| simulations predicting that the W+ and W~ 
bosons are both predominantly left-handed in W+jets events at high pj. An NLO QCD calcu- 
lation [36J has demonstrated that the predicted polarization fractions are stable with respect to 
QCD corrections. As discussed in later sections, this detailed knowledge of the W-boson polar- 
ization provides key information for measuring the SM backgrounds using control samples in 
data. 

4 Event preselection 

Table [I] summarizes the main variables and requirements used in the event preselection, which 
is designed to be simple and robust. Except where noted, a common set of preselection require- 
ments is used by each of the three analysis methods. Events are required to have at least one 
good reconstructed primary vertex, at least three jets (Lp method and ANN method) or four 
jets (LS method), and exactly one isolated muon or exactly one isolated electron. These basic 
requirements select an event sample that is dominated by genuine, single-lepton events from 
SM processes. 

The primary vertex must satisfy a set of quality requirements, including |zpy| < 24 cm and 
ppv < 2 cm, where zpy and ppy are the longitudinal and transverse distances of the primary 
vertex with respect to the nominal interaction point in the CMS detector. 

Jets are reconstructed offline using the anti-Zcp clustering algorithm fl37) with a distance param- 
eter of 0.5. The particle four- vectors reconstructed by the CMS particle-flow algorithm p3ll24| , 
are used as inputs to the jet clustering algorithm. The particle-flow algorithm combines in- 
formation from all CMS sub-detectors to provide a complete list of long-lived particles in the 
event. Corrections based on simulation are applied to the jet energies to establish a uniform 
response across the detector and a first approximation to the absolute energy scale ||38|| . Addi- 
tional jet energy corrections are applied to the data using measurements of energy balance in 
dijet and photon + jet control samples in data. These additional corrections take into account 
residual differences between the jet energy scale in data and simulation. The effect of pileup 
was significant during much of the data-taking period. Extra energy clustered into jets due 
to pileup is taken into account with an event-by-event correction to the jet momentum four- 
vectors. Jet candidates are required to satisfy quality criteria that suppress noise and spurious 
energy deposits in the calorimeters. The performance of jet reconstruction and the corrections 
are described in Ref. |38j. In this analysis, reconstructed jets are required to satisfy px > 40 GeV 
and \rj \ < 2.4. The £x vector is defined as the negative of the vector sum of the transverse 
momenta of all the particles reconstructed and identified by the particle-flow algorithm. 

In the muon channel, the preselection requires a single muon candidate II39II satisfying px(p) > 
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Table 1: Main preselection requirements. The term lepton designates either an electron or a 
muon. Definitions of the quantities and further details are given in the text. 
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Jet t] range 


|f/| < 2.4 


Number of jets 


> 3 (L P and ANN methods), 




> 4 (LS method) 


Lepton px threshold 


> 20GeV 


Muon rj range 


\V\ < 2- 1 


Muon isolation (relative) 


< 0.10 


Electron rj range 


7/| < 1.442, 1.56 < 7/1 < 2.4 


Electron isolation (relative) 


< 0.07 (barrel) 




< 0.06 (endcaps) 


Lepton pj thresh, for veto 


> 15GeV 



20GeV and \tj\ < 2.1. Several requirements are imposed on the elements that form the muon 
candidate. The reconstructed track must satisfy quality criteria related to the number of hits in 
the pixel, strip, and muon detectors, and it must have an impact parameter do in the transverse 
plane with respect to the beam spot satisfying \do\ < 0.02 cm and an impact parameter d z with 
respect to the primary vertex along the z direction satisfying \d z \ < 1.0 cm. 

To suppress background in which the muon originates from a semileptonic decay of a hadron 
containing a bottom or charm quark, we require that the muon candidate be spatially isolated 
from other energy in the event. A cone of size AR = 0.3 is constructed around the initial 
muon momentum direction in tj-(p space. The muon combined isolation variable, J comb = 
Ear<o.3(^t + Pt)/ is defined as the sum of the transverse energy Ej (as measured in the elec- 
tromagnetic and hadron calorimeters) and the transverse momentum pj (as measured in the 
silicon tracker) of all reconstructed objects within this cone, excluding the muon. This quan- 
tity is used to compute the combined isolation relative to the muon transverse momentum, 
jcomb = JC omb /pT (^ which ig require( j to sa tisfy 7 r c e ° mb < 0.1. 

Electron candidates [40 1 are reconstructed by matching energy clusters in the ECAL with tracks 
in the silicon tracking system. Candidates must satisfy pj > 20 GeV and \t]\ < 2.4, excluding 
the barrel-endcap transition region (1.442 < \tj\ < 1.56). Quality and photon-conversion rejec- 
tion requirements are also imposed. The relative isolation variable, defined in a manner similar 
to that in the muon channel, must satisfy JJr°j mb < 0.07 in the barrel region and 7 r c e ° 1 mb < 0.06 in 
the endcaps. The requirements on do and d z are the same as those used in the muon channel. 

The preselection requirements have a large effect on the sample composition. The lepton isola- 
tion requirement is critical for the rejection of QCD multijet processes, which have very large 
cross sections but are reduced to a low level by the isolation and the other preselection require- 
ments. While many lepton candidates are produced in the semileptonic decays of hadrons 
containing b or c quarks, from n and K decays in flight, and from misidentification of hadrons, 
the vast majority of these candidates are either within or near hadronic jets. The background 
from W+jets events (primarily from W — > ev or W — > y.v, but also W — > tv) is initially also 
very large. This contribution is heavily suppressed by the three- or four-jet requirement. De- 
pending on the particular signal region, either tt or W+jets production emerges as the largest 
contribution to the background in the sample of events with moderate to large lj. 
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5 Lepton Spectrum method 



Events with a second isolated-lepton candidate satisfying the criteria listed in Table [T] are ve- 
toed. This requirement not only suppresses SM background, but also minimizes the statistical 
overlap between the event sample used in this search and those used in multilepton searches. 
However, tt events with dileptons can still be present, and this contribution must be deter- 
mined, particularly because the presence of two neutrinos in the decay chains can result in 
large values of The background involving W — > tv decays, both from tt events and from 
direct W production, must also be determined. To help suppress the dilepton background, the 
requirements on the veto leptons are somewhat looser than those on the signal lepton. For both 
muons and electrons, the pj threshold is pj > 15 GeV, the isolation requirement is I^ mb < 0.15, 
and the impact parameter requirement is \do\ <0.1 cm (the d z requirement is kept the same as 
for the signal lepton). In addition, some of the quality requirements for both the muon and 
electron are loosened. 

Further event selection requirements are used in the individual background estimation meth- 
ods described in Sections |5}|6} and[7| The methods use the quantity Hj, which is defined as the 
scalar sum of the transverse momenta of particle-flow jets j with p J T > 40 GeV and < 2.4, 

H T = X> ; T . (2) 
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The three background determination methods presented in the following three sections use 
different approaches to estimating the SM backgrounds using control samples in data. In Sec- 
tion [9j we compare the results of the different methods and make some observations about 
their features. 

5 Lepton Spectrum method 

5.1 Overview of the Lepton Spectrum method 

This section describes the Lepton Spectrum (LS) method, which is named for the technique 
used to determine the dominant background source: genuine, single-lepton processes. Such 
processes account for about 75% of the total SM background in the signal regions and arise 
primarily from tt, single-top, and W+jets events. Their contribution to the #j distribution is 
estimated by exploiting the fact that, when the lepton is produced in W-boson decay, the #j 
distribution is fundamentally related to the lepton pj spectrum, unlike the £t for many SUSY 
models. A more detailed description of the Lepton Spectrum method is given in the refer- 
ences bush. 

Non-single-lepton backgrounds are also determined using control samples in the data. Such 
events arise mainly from (i) tt dilepton events, in which zero, one, or both of the leptons is a 
t and (ii) tt and W+jets events with a single x — > (ji, e) decay. Background from QCD mul- 
tijet events is expected from simulation to be very small. However, the uncertainties in such 
simulations are difficult to quantify, because the QCD multijet background in the phase space 
relevant to this analysis arises from extreme tails of processes with very large cross sections. 
We therefore use control samples in data to measure the QCD multijet background. Simulated 
event samples are used for the determination of the Z+jets contribution, which is estimated 
with sufficient precision to be below one event for most of the signal regions. 

The signal regions are defined with three thresholds in Hj (Hj > 500 GeV, Hj > 750 GeV, and 
Hj > 1000 GeV) and four bins in l T (250 < JS T < 350 GeV, 350 < £ T < 450 GeV, 450 < j£ T < 
550 GeV, and £ T > 550 GeV). 



5.2 Estimation of single-lepton backgrounds 
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Figure 1: Distributions of muon pj vs. £x in the p channel for simulated tt and W+jets events 
(left) and for the LM6 SUSY benchmark model (right). 

5.2 Estimation of single-lepton backgrounds 

The physical foundation of the Lepton Spectrum method is that, when the lepton and neutrino 
are produced together in two-body W decay (either in tt or in W+jets events), the lepton pj 
spectrum is directly related to the £x spectrum. The lepton and the neutrino share a common 
Lorentz boost from the W rest frame to the laboratory frame. As a consequence, the lepton spec- 
trum reflects the pj distribution of the W, regardless of whether the lepton was produced in a 
top-quark decay or in a W+jets event. With suitable corrections, discussed below, the lepton pj 
spectrum can therefore be used to predict the %j spectrum for SM single-lepton backgrounds. 

The £j distribution in many SUSY models is dominated by the presence of two LSPs. In con- 
trast to the SM backgrounds, the £x and lepton pj distributions in SUSY processes are therefore 
nearly decoupled. The £x distribution for such models extends to far higher values than the 
lepton spectrum. Figure [T] shows the relationship between the lepton-px and £x distributions 
in the laboratory frame for two simulated event samples: (i) the predicted SM mixture of tt and 
W+jets events and (ii) the SUSY LM6 benchmark model. When taken from data, the upper-left 
region in Fig. [I] (left) provides the key control sample of high-px leptons from SM processes. 
This region typically has very little contamination from SUSY events, which populate the high- 
£x region but have relatively low lepton pj values. 

The lepton pj spectrum is measured with a muon control sample defined by the preselection 
criteria and the Hj requirements. Unlike the signal region, no £x requirement is applied, be- 
cause even a modest one (£x > 25 GeV) would bias the high end of the lepton pj spectrum, 
which is critical for making the background prediction. Only muon events are used as a con- 
trol sample, because the QCD multijet background is significant in the low-^x region of the 
electron sample. The number of events that are common to both the control sample and the 
signal region is small. For example, the overlap as measured in simulated tt events is 3.6% for 
Hj > 750 GeV, #x > 250 GeV, and pj > 250 GeV. Because no £x requirement is placed on 
the muon control sample, a small amount of QCD background remains and must be measured 
and subtracted. The scaling from the muon to the electron samples is obtained by fitting their 
ratio in the data over the range 60 < < 250 GeV, with systematic uncertainties evaluated by 
varying the fit range. The resulting correction factor is N(e)/N(p.) = 0.88 ± 0.03 ± 0.03, where 
the uncertainties are statistical and systematic, respectively. 
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5 Lepton Spectrum method 



To use the lepton spectrum to predict the Ej spectrum in single-lepton SM background pro- 
cesses, three main issues must be understood: (i) the effect of the W-boson polarization in both 
tf and W+jets events, (ii) the effect of the applied lepton pj threshold, and (iii) the difference 
between the experimental resolutions on the measurements of lepton pj and Ej. 

The status of theoretical and experimental knowledge of W-boson polarization in tf and in 
W+jets events is discussed in Section|3j The helicity zero polarization state results in a forward- 
backward symmetric angular distribution of the lepton and the neutrino in the W rest frame 
(with respect to the W momentum direction), leading to identical lepton and neutrino spectra in 
the laboratory frame. In contrast, the helicity ±1 states result in angular asymmetries that lead 
to somewhat different lepton and neutrino pj spectra in the laboratory frame. These effects are 
taken into account by applying correction factors obtained from simulation to the measured 
lepton spectrum, with uncertainties as described in Section [HJ 

The second key issue in the Lepton Spectrum method is the effect of the threshold (pj > 
20 GeV) applied to the leptons in both the signal and control samples. Because of the anti- 
correlation between the lepton pj and the Ej arising from non-zero W-boson helicity states, the 
threshold requirement removes SM background events in the high-^j signal region but not the 
events in the control sample with high-px muons that are used to predict the high tail of the 
Ej spectrum. For the tt background, this effect partially compensates for the bias from the W 
polarization. For W+jets events, in contrast, the polarization effects for W + and W~ approxi- 
mately cancel, but the lepton px threshold shifts the predicted yield upward. Correction factors 
from simulation are used to account for these effects (as well as for polarization effects), which 
are well defined and understood. 

Finally, the resolution on the reconstructed Ej is poorer than that for the lepton px, so the Ej 
spectrum is somewhat broadened with respect to the prediction from the lepton spectrum. We 
measure Ej resolution functions in the data using QCD multijet events obtained with a set of 
single-jet triggers spanning the range from Ej > 30 GeV to Ej > 370 GeV. These resolution 
functions, or templates, quantify the Ej resolution as a function of the number of jets and the 
Hj of the event. These templates are used to smear the measured lepton momenta. Because 
the templates are taken from data, they include not only the intrinsic detector resolutions, but 
also acceptance effects. The overall effect of the smearing is modest, changing the background 
prediction by 5-15%, depending on the Ej threshold applied. 

The raw background predictions for the single-lepton background are corrected to account for 
the effects described above, as well as for the small contamination of the single-lepton control 
sample arising from dilepton and single-T events with high-px leptons. These backgrounds are 
measured separately, as described below. The overall correction factor is defined such that the 
single-lepton prediction in a given signal region in simulation matches the yield from single- 
lepton processes. 

The predicted single-lepton background yield varies from about 150 events for the signal region 
with 250 < Ej < 350 GeV and Hj > 500 GeV to about 3 events for the region with E T > 
550 GeV and Hj > 1000 GeV. These predictions, as well as the expectations from simulation, 
are presented in Tables|2j|3j and|4]and discussed in more detail in Section 5.4 



5.3 Estimation of non-single-lepton backgrounds 

The non-single-lepton backgrounds include dilepton events in several categories, events with 
W — > tv followed by r — > £ decays (in either tf or W+jets events), and QCD multijet processes. 
These subdominant backgrounds are estimated using control samples in data, in conjunction 



5.3 Estimation of non-single-lepton backgrounds 
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Table 2: Event yields for the Lepton Spectrum method for Hj > 500 GeV. The upper part of 
the table gives the background predictions that are based on simulated (MC) event samples 
and the yield for the SUSY signal points LM3 and LM6. The lower part gives the backgrounds 
predicted using control samples in the data (data-driven prediction). The actual yield observed 
in data is given at the bottom, with the separate muon and electron yields given in parentheses 
(Nji/Ne) after the total yield. The uncertainties on the background predictions are statistical 
and systematic. The MC yields are not used in setting limits and are included only for reference. 
The uncertainties on the MC yields are statistical only. 



#t range [GeV] 


[250, 350) 


[350, 450) 


[450, 550) 


>550 






MC yields 






1 1 


146.7 ± 2.1 


34.8 ± 1.1 


8.5 ± 0.6 


2.9 ± 0.3 


Dilepton 


19.9 ± 0.5 


3.8 ± 0.2 


0.7 ± 0.1 


0.3 ±0.1 


1 T 


30.6 ± 0.9 


7.9 ± 0.5 


2.1 ± 0.3 


0.8 ± 0.2 


Z+jets 


1.3 ±0.8 


<0.1 


<0.1 


<0.1 


Total (MC) 


198.6 ± 2.5 


46.5 ± 1.2 


11.3 ±0.6 


4.0 ± 0.4 


SUSY LM3 (MC) 


266.3 ± 3.7 


91.0 ± 2.2 


23.3 ± 1.1 


9.9 ± 0.7 


SUSY LM6 (MC) 


23.4 ± 0.3 


20.0 ±0.3 


13.4 ±0.2 


10.8 ±0.2 




Data-driven prediction 






It 


109 ± 13 ± 18 


32.0 ± 7.5 ± 5.8 


3.9 ± 2.7 ± 1.2 


3.1 ± 2.3 ± 1.0 


Dilepton 


15.8 ±1.9 ±1.8 


3.0 ± 0.9 ± 0.5 


0.5 ± 0.3 ± 0.2 


0.1 ± 0.2 ± 0.2 


1 T 


33.0 ± 1.8 ± 1.7 


8.9 ± 1.0 ± 0.5 


2.1 ±0.5 ±0.2 


1.1 ± 0.3 ± 0.2 


QCD 


0.0 ± 1.0 ± 1.2 


0.0 ± 1.0 ± 1.2 


0.0 ± 1.0 ± 1.2 


0.0 ± 1.0 ± 1.2 


Z+jets (MC) 


1.3 ± 0.8 ± 1.3 


<0.1 


<0.1 


<0.1 


Total (predicted) 


159 ± 14 ± 18 


44.0 ± 7.7 ± 6.0 


6.6 ±2.9 ±1.7 


4.3 ± 2.6 ± 1.6 


Data (observed) 


163 (84, 79) 


46 (21, 25) 


9 (8, 1) 


2(1,1) 



Table 3: Event yields for the Lepton Spectrum method for Hj > 750 GeV. Further details are 
given in the Table [2] caption. 



$T range [GeV] 


[250, 350) 


[350, 450) 


[450, 550) 


>550 






MC yield 






It 


47.3 ± 1.2 


14.9 ± 0.7 


5.4 ± 0.4 


2.7 ±0.3 


Dilepton 


8.2 ± 0.4 


2.3 ± 0.2 


0.6 ± 0.1 


0.3 ±0.1 


1 T 


9.2 ± 0.5 


3.0 ± 0.3 


1.2 ±0.2 


0.7 ±0.2 


Z+jets 


0.7 ±0.6 


<0.1 


<0.1 


<0.1 


Total (MC) 


65.4 ± 1.5 


20.2 ± 0.8 


7.2 ± 0.5 


3.6 ± 0.4 


SUSY LM3 (MC) 


114.6 ± 2.5 


47.1 ± 1.6 


16.1 ±0.9 


8.6 ± 0.7 


SUSY LM6 (MC) 


14.9 ± 0.3 


13.8 ±0.2 


10.3 ±0.2 


9.8 ± 0.2 




Data-driven prediction 






1 i 


41.7 ±8.7 ±5.4 


11.7 ± 5.0 ± 1.9 


2.6 ± 2.3 ± 0.6 


3.1 ± 2.4 ± 0.8 


Dilepton 


5.9 ± 1.1 ± 0.7 


1.3 ± 0.5 ± 0.2 


0.5 ±0.2 ±0.1 


0.1 ± 0.1 ± 0.3 


1 T 


9.6 ± 0.9 ± 0.6 


3.1 ± 0.6 ± 0.3 


1.1 ±0.3 ±0.2 


0.8 ±0.2 ±0.1 


QCD 


0.0 ± 0.2 ± 0.4 


0.0 ± 0.2 ± 0.4 


0.0 ± 0.2 ± 0.4 


0.0 ± 0.2 ± 0.4 


Z+jets (MC) 


0.7 ±0.6 ±0.7 


<0.1 


<0.1 


<0.1 


Total (predicted) 


57.9 ± 8.9 ± 5.6 


16.2 ± 5.0 ± 2.0 


4.2 ± 2.4 ± 0.8 


4.0 ± 2.4 ± 1.0 


Data (observed) 


48 (27, 21) 


16 (7, 9) 


5 (4, 1) 


2(1,1) 
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Table 4: Event yields for the Lepton Spectrum method for Hj > 1000 GeV. Further details are 
given in the Table [2] caption. 



J?T range [GeV] 


[250, 350) 


[350, 450) 


[450, 550) 


>550 






MC yield 






11 


13.4 ± 0.6 


4.8 ± 0.4 


2.1 ± 0.3 


1.3 ±0.2 


Dilepton 


2.7 ± 0.2 


1.0 ±0.1 


0.3 ± 0.1 


0.2 ±0.1 


1 T 


2.1 ± 0.2 


0.7 ±0.1 


0.5 ± 0.1 


0.4 ±0.1 


Z+jets 


0.5 ± 0.5 


<0.1 


<0.1 


< 0.1 


Total (MC) 


18.8 ± 0.9 


6.4 ± 0.5 


2.9 ± 0.3 


1.9 ± 0.2 


SUSY LM3 (MC) 


38.1 ± 1.4 


18.3 ± 1.0 


7.0 ± 0.6 


5.5 ± 0.5 


SUSY LM6 (MC) 


7.0 ± 0.2 


6.0 ± 0.2 


4.6 ±0.1 


5.2 ± 0.2 




Data-driven prediction 
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11.7 ±4.6 ± 1.8 


5.5 ± 3.6 ± 1.0 


2.0 ± 2.2 ± 0.6 


3.1 ± 2.3 ± 1.0 


Dilepton 


1.2 ±0.6 ±0.1 


0.4 ± 0.4 ± 0.1 


0.2 ±0.2 ±0.1 


0.1 ±0.2 ±0.2 


1 T 


3.0 ± 0.5 ± 0.5 


0.9 ± 0.3 ± 0.2 


0.4 ± 0.2 ± 0.2 


0.8 ± 0.2 ± 0.2 


QCD 


0.0 ±0.1 ±0.1 


0.0 ± 0.1 ± 0.1 


0.0 ±0.1 ±0.1 


0.0 ±0.1 ±0.1 


Z+jets (MC) 


0.5 ± 0.5 ± 0.5 


<0.1 


<0.1 


< 0.1 


Total (predicted) 


16.4 ± 4.7 ± 1.9 


6.8 ± 3.6 ± 1.0 


2.6 ± 2.2 ± 0.6 


4.0 ± 2.4 ± 1.0 


Data (observed) 


14 (7, 7) 


4(1,3) 


(0, 0) 


2(1,1) 



with information from simulation. The contribution from Drell-Yan and Z+jets is very small 
and is estimated directly from simulation. 

Dilepton background events (including the t as one of the leptons) contain at least two neutri- 
nos, so these events can be important in the tails of the £x distributions. These backgrounds are 
divided into the following categories: (i) 2.1 events with one lost or ignored lepton (£ = e, 
(ii) I + T events with x — > hadrons, and (iii) £ + T events with X — > lepton. A lost lepton is one 
that is either not reconstructed or is out of the detector acceptance. An ignored lepton is one 
that is reconstructed but fails either the lepton-identification requirements (including isolation) 
or the px threshold requirement. 

To estimate the background from dilepton events with lost or ignored leptons, we compute 
the ratio of the combined yield of dilepton events in the ee, efi, and jiji channels in data to 
the corresponding combined yield in simulated event samples. This ratio, which is 0.91 ± 0.07 
for Hj > 500 GeV, 0.93 ± 0.15 for Hj > 750 GeV, and 0.87 ± 0.37 for Hj > 1000 GeV, is used 
to rescale the £x distribution of dilepton events that appear in the signal region in simulation. 
(Events within 20 GeV of the nominal Z mass are excluded in the e + e and channels.) 
This approach is used because the dilepton control sample in data is small, and using it to 
obtain the shapes of £x distributions would result in large statistical uncertainties. For all $j 
bins above 250 GeV, the predicted yield from this background contribution is less than 6 events, 
and for all bins above 350 GeV, the yield is at or below 1 event. The £j distribution associated 
with the reconstructed dilepton events in data is well described by the simulation. 

Dilepton events can also involve x decays, either x — > hadrons or x — > i. The Ifcj distributions 
in the dilepton events in data, when suitably modified to reflect the presence of a leptonic or 
hadronic x decay, provide an accurate description of the shape of the £t distribution of these 
backgrounds. Thus, to estimate the shape from the x — > hadrons background, we effectively 
replace a lepton in a reconstructed dilepton event with a hadronic X jet. Both hadronic and 
leptonic x response functions are used, providing a probability distribution for a x to produce 
a jet or a lepton with a given fraction px(jet)/ pj(x) or pj(£) /pj(x). These response functions, 



5.3 Estimation of non-single-lepton backgrounds 
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Figure 2: Predictions for dilepton and T — > £ backgrounds after requiring Hj > 750 GeV: con- 
trol samples in data (green points with error bars) vs. MC predictions (black solid histogram) 
for (left) dilepton background and (right) T — > £ background. The MC prediction has been 
scaled to the integral of the data prediction. 

obtained from simulation, are computed in bins of Pt(t). This procedure can change the total 
number of jets above threshold in the event, as well as other properties such as Hj and £t, 
which are recalculated. Simulated event samples are used to determine, for these processes, 
the ratio r, = N' ieed / N contro i of the number of events observed in the single-lepton channel to 
the number of events in the control sample, as a function of Ifcj. This procedure effectively 
normalizes all such contributions to the control samples in data. For all £j bins above 250 GeV, 
the number of dilepton events with a t — > hadrons decay is predicted to be about 7 events or 
less and is much smaller in the higher £x bins. The number of dilepton events with a T — > £ 
decay is predicted to be less than 3 events for all bins above 250 GeV and is much smaller in 
the higher £x bins. 

Estimates for the T — > £ single-lepton backgrounds from tt and W+jets processes are based 
on a procedure similar to that used for the dilepton backgrounds, but in this case the single- 
lepton sample itself is used as the control sample. The £j distribution obtained by applying the 
t — > £ response function to the data is rescaled by a ratio from simulation that gives the yield 
of t — > £ background events divided by the yield of events in the single-lepton control sample, 
as a function of Ifcj. The number of background events from the single T — > £ contribution 
falls from 33 for Hj > 500 GeV and 250 < £ T < 350 GeV to 1.1 event for Hj > 500 GeV and 
It > 550 GeV. 

The background predictions in data are shown in Fig. |2j where the expectation based on sim- 
ulation is shown for comparison. The total predicted dilepton plus single T — > £ background 
yield ranges from about 50 events for Hj > 500 GeV and 250 < ?t < 350 GeV to about 1 event 
for Hj > 1 TeV and > 550 GeV. All of these predictions, as well as the expecations from sim- 
ulation, are presented in Tables |2j|3j and|3J which are discussed in more detail in Section 5.4 



Background from QCD multijet events is suppressed to a level well below 1 event in nearly all 
signal regions, as shown in Tables |2j|3j and|4] The QCD multijet background is determined by 
first defining a control sample with small missing transverse momentum (^x < 50 GeV) and 
with a lepton impact parameter relative to the beam spot |c?o| > 0.02 cm. These requirements 



12 



6 Lepton Projection method 



select a sample with little contamination from other SM processes such as tt and W+jets pro- 
cesses. Using this control sample, we measure the shape of the distribution in the combined 
relative isolation variable, /^°j mb (see Section 5j. The shape of this distribution has very little 
correlation to #j or to the lepton impact parameter (do), and so can be applied in the high-^x 
signal regions. For each signal region in the data, we determine the background at low values 
of J^° 1 mb by first scaling the measured QCD multijet background shape in the relative isolation 
variable to the high-7™ 1 mb sideband of the signal region. The shape is then used to extrapolate 
the yield to the low- ^° 1 mb signal region. In the high-^x signal regions, some non-QCD SM 
background can be present at high I^j mb / where the QCD background shape is normalized. We 
therefore subtract the estimated background from tt, W+jets, and Z+jets from this region. These 
yields are taken from simulation, with systematic uncertainties determined from a comparison 
with a control region in the data. 

5.4 Results from the Lepton Spectrum method 

Tables |2j |3} and [4] compare the background yields predicted from the control samples in data 
with the yields obtained directly from simulation for Hj > 500 GeV, Hj > 750 GeV, and 
Hj > 1000 GeV, respectively. We observe that the single-lepton background is the dominant 
contribution in all regions. The various sources of uncertainties associated with these back- 
ground determinations are discussed in Section [9] Finally, the yields observed in the signal 
regions in the data, which are listed at the bottom of each table, are consistent with the total 
background predictions based on the control samples. Thus, we observe no evidence for any 
excess of events in the data above the SM contributions. 

Figure [3] shows the distributions in data for the combined muon and electron channels, 
with all of the selection requirements, except that on £j itself. The distributions are shown 
for Hj > 500 GeV, Hj > 750 GeV, and Hj > 1000 GeV, on both linear and logarithmic scales. 
The predicted £x distribution (green-bar histogram) is a sum over three sources: single-lepton 
backgrounds (from tf, single-top, and W+jets events), dilepton background from tf, and single- 
t events (from both tf and W+jets processes). The vertical span of the green bar corresponds 
to the statistical uncertainty on the background prediction. (The systematic uncertainties are 
computed in wider bins used for setting the limits and are given in Tables |2j |3j and |4]) In 
each signal region, the blue histogram shows the contribution from the dilepton and single-r 
backgrounds only. It is evident that the single-lepton background is dominant in all cases. The 
£x distributions for the SUSY benchmark models LM3 and LM6 are overlaid (not summed) for 
comparison. Systematic uncertainties and the interpretation are presented in Section [9] 

6 Lepton Projection method 

6.1 Overview of the Lepton Projection method 

The Lepton-Projection (LP) method uses the difference between SM and SUSY processes in the 
correlation of the lepton transverse momentum and the missing transverse momentum. As 
previously discussed, in the SM processes the corresponds to the neutrino in the decay of 
the W boson, either in W+jets or in tt events. The kinematics of W decays are dictated by the 
V— A nature of the W coupling to leptons and the helicity of the W boson, as discussed in 
Section[3] Since W bosons that are produced with high transverse momentum in W+jets events 
exhibit a sizable left-handed polarization, there is a significant asymmetry in the pj spectra of 
the neutrino and charged lepton. A smaller asymmetry is expected in W bosons from t quark 
(t antiquark) decays, which yield W bosons which are predominantly longitudinally polarized 
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Figure 3: Lepton Spectrum method: observed %j distributions in data (filled points with er- 
ror bars) compared with predicted £x distributions (green bars) in the combined electron and 
muon channels, on linear (left) and logarithmic (right) scales. Three different Hj thresholds are 
applied: H T > 500 GeV (upper row), H T > 750 GeV (middle row), and H T > 1000 GeV (lower 
row). 



14 



6 Lepton Projection method 



with smaller left-handed (right-handed) components for W + (W ). 



We have measured the fraction of the helicity states of the W boson using an angular analysis 
of leptonic W decays l|2l) . Since the total momentum of the W boson in these decays, and 
therefore its center-of-mass frame, cannot be accurately determined because the momentum 
of the neutrino along the beam axis cannot be measured, an observable that depends only on 
transverse quantities is used. A variable that is highly correlated with the cosine of the polar 
angle in the center-of-mass frame of the W boson is the "lepton projection variable": 



p T (l) • Pt(W) 
IPt(W)| 2 ' 



(3) 



where pi(£) is the transverse momentum of the charged lepton and px(W) is the transverse 
momentum of the W boson. The latter quantity is obtained from the vector sum of the electron 
transverse momentum and the missing transverse momentum in the event. 

Since SUSY decay chains result in large values of Ej, and often result in relatively low values 
of the lepton momentum as well, the Lp distribution for SUSY events tends to peak near zero, 
whereas W+jets and tt yield a broad range of Lp values. This behavior is illustrated in Fig. |Ij 
which compares the Lp distribution from both SM processes and from two representative SUSY 
benchmark points (LM3 and LM6). 
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Figure 4: Distribution of Lp in SUSY and standard model processes from simulation. Left: all 
distributions are normalized to the integrated luminosity. The different contributions from SM 
processes are shown, whereas for SUSY two benchmark points, LM3 and LM6, are displayed. 
Right: the same distributions normalized to unity. The SM distribution is the sum of all the 
individual SM processes shown in the left pane. The quantity Sj P = pj(() + is discussed in 
the text. 



In the Lp method, two regions in Lp are defined: the region with Lp < 0.15 is used as the signal 
region; the region with Lp > 0.3 is used as the control region, i.e., a sample that is depleted 
in the signal expected and is instead dominated by SM processes. These regions are selected 
using simulated event samples of W+jets, Z+jets, and tt, that are collectively referred to as 
electroweak (EWK) processes in what follows, as well as with simulated SUSY events with 
SUSY particle masses near the region currently under exploration. 



6.2 Background estimation in the L P method 
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6.2 Background estimation in the L P method 

The key ingredient of the analysis is the estimate of the number of events in the signal region 
from the SM processes. We define a translation factor, 

N MC (L P < 0.15) 
KCS N MC (L P >0.3)' W 

which is the ratio of the number of events in the signal and control regions for the EWK pro- 
cesses. The translation factor is obtained from MC simulation of the EWK processes, and the 
uncertainties on this factor are included in the systematic uncertainty of the background esti- 
mate. In the case of muons, where the background from QCD multijets is negligible, the total 

number of events predicted from SM processes in the signal region, N P j^ d (Lp < 0.15), can 
be determined directly from the number of events observed in the data in the control region, 
Ndata^P > 0.3): 

N s p ^ ed (Lp < 0.15) = R CS • N data (L P > 0.3). (5) 

In the case of the electrons, the presence of events from QCD multijet processes necessitates an 
independent evaluation of this background prior to the application of the translation factor for 
EWK processes. 

The number of events estimated with this method is then compared to the number of events 
observed in the data in the signal region, N^ ata (Lp < 0.15), for indications of an excess of 
events over the SM expectation. The analysis is performed in different regions of the event 
mass scale. To characterize the latter without affecting the correlation of the charged lepton 
and the neutrino in SM events, the scalar sum of the lepton transverse momentum and the 
missing transverse momentum, Sj P , is used: Sj P = pj(£) + tj- For W decays, Sj P ~ pr(W) 
at large values of px(W). 

In order to make the search optimization less dependent on the unknown energy scale of a 

new physics signal, the analysis is performed in disjoint ranges of Sj P and the results in these 
ranges are combined. In addition, the selection is also binned in a second dimension, the Hj 
variable, defined in Eq. 

As indicated in Table [TJ the event selection used in this analysis is slightly different from the 
corresponding one in the LS analysis. To increase the sensitivity to SUSY decays, this analysis 
requires three or more jets. While this results in a significant increase in W+jets events, the 
additional SM background is mostly concentrated in the control region in Lp. 

The event yields in the muon and electron channels, as predicted from simulation, are shown 
in Table [5] As discussed previously, the dominant backgrounds to the lepton plus jets and ?x 
signature arise from the production and decay of W+jets and tt. The production of single W 
bosons in association with jets, and with large transverse momenta, is in general the larger of 
the two, especially at lower jet multiplicities. The majority of the tt background arises from 
semi-leptonic tt decays, with fully leptonic tt decays in which a lepton is either ignored or not 
reconstructed contributing about 20% of the total tt background. 

A source of background, which is not listed in Table |5j stems from QCD multijet events in 
which a jet is misreconstructed as a lepton. The simulation indicates that the magnitude of this 
background is small in the control region and negligible in the signal region. Nevertheless, 
since the uncertainties in simulating these backgrounds can be significant, we use control data 
samples to estimate the background in the muon and electron channels. 
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Table 5: Expected event yields in the signal region (Lp < 0.15) from simulation. These yields 
are for Hj > 500 GeV. These MC values are only listed for illustration purposes. 



L P < 0.15 


Muons: Sj P range [GeV] 


Electrons: Sj P range [GeV] 




[250-350] 


[350-450] 


[450-oo] 


[250-350] 


[350-450] 


[450-oo] 


tt(£) 


50.0 ±1.0 


15.3 ±0.5 


4.8 ±0.3 


37.9 ± 0.8 


11.0 ±0.4 


3.6 ±0.2 


tt (££) 


12.4 ±0.4 


3.9 ±0.2 


1.2 ±0.1 


10.4 ±0.4 


2.9 ±0.2 


0.8 ±0.1 


W 


66.2 ±2.0 


35.6 ±1.4 


26.0 ±1.2 


48.9 ±1.7 


24.2 ±1.2 


20.9 ±1.1 


Z 


2.1 ±1.0 


0.4 ± 0.4 


0.0 ±0.2 


1.4 ±0.8 


0.0 ±0.2 


0.0 ±0.2 


Total MC 


130.8 ±2.4 


55.3 ±1.6 


32.0 ±1.3 


98.6 ±2.1 


38.1 ±1.3 


25.3 ±1.1 


LM3 


136.8 ±3.8 


89.1 ±3.1 


53.9 ± 2.4 


111.7 ±3.4 


70.8 ± 2.7 


47.0 ± 2.2 


LM6 


8.4 ±0.2 


11.0 ±0.2 


24.9 ± 0.3 


6.7 ±0.2 


8.5 ±0.2 


20.5 ±0.3 



To estimate the background from QCD multijets in the muon final state, we use the relative 
combined isolation, IJ^jL,/ of the muon. Multijet events are expected to populate the region 
at high values of I*omb' whereas muons from SUSY decays are isolated and thus have low 
values of I^mb - ^ e em pl°y 311 additional control data sample, which is specially selected to be 
enriched in QCD multijets, to determine the ratio of multijet events at low values of the relative 
isolation. Using this ratio and the number of multijet events expected in the control region of 
the sample passing the preselection requirements, we estimate the background from multijet 
events in the signal region to be always smaller than 1% of the EWK backgrounds. This level 
of background is negligible and is thus ignored in what follows. 

The main sources of electrons in QCD multijet events are misidentified jets and photon con- 
versions. This background is expected to be more substantial than the corresponding one in 
the muon sample, and its estimate exhibits a large dependence on the details of the simulation. 
For this reason, we estimate this background from control samples in data. The method relies 
on the inversion of one or more of the electron identification requirements in order to obtain a 
sample of anti-selected events, which is dominated by jets misidentified as electrons. We find 
that the inversion of the requirements on the spatial matching of the calorimeter cluster and 
the charged-particle track in pseudorapidity and azimuth leaves the relative fraction of the dif- 
ferent background sources in QCD multijets unchanged. Moreover, to increase the number of 
events in this control sample, the requirements on do and d z are removed, while the isolation 
requirement is loosened. These changes to the event selection have a negligible effect on the 
Lp distribution in the data. In the simulated event samples, it is found that the Lp distribution 
from the control sample events provides a good description of the corresponding distribution 
from QCD background passing all selection requirements. 

The Lp distribution obtained with this control sample is used as a template to fit, along with 
the Lp distribution from EWK processes, the Lp distribution in the data. In this fit, the EWK 
template is taken from simulation. This approach, which provides a template obtained from 
data for the QCD contamination, was applied in the measurement of the polarization of high- 
pj W bosons [21]. The fit is performed in the control region (Lp > 0.3), where the possible 
presence of signal is highly suppressed. The numbers of QCD and EWK events obtained by 
the fit are used to estimate the total SM contamination in the signal region (Lp < 0.15). The 
method for estimating the number of SM events expected in the signal region is applied in 
each range of Sj P and Hj. 

The method for estimating the SM expectation in the signal region is checked using two dif- 



6.3 Results of the L P method 
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Figure 5: Fit results on data for 150 < < 250 GeV, in the muon (left) and electron (right) 
search samples. The fit is performed in the control region (Lp > 0.3) and the result is extrapo- 
lated into the signal region (Lp < 0.15). 



ferent control samples, where both the fit and signal regions have a negligible expected SUSY 
yield. The first sample is defined as all events satisfying the preselection requirements but con- 
fined to low values of Sj P : 150 < Sj P < 250 GeV. The method described above is employed 
to predict the number of events expected in the signal region for both muons and electrons. 
This prediction is found to be fully consistent with the number of events observed in the data 
in signal region. The results of the fits and the yields of QCD and EWK events in the region 

of low S j p ( < 250 GeV) are displayed in Fig. [5] for the electron and muon samples. As can be 
seen in Fig. [5} the QCD contamination in the signal region, Lp < 0.15, is negligible, as expected, 
since low values of Lp favor events with low-px leptons and high J£j- The second sample, used 
only for events with muons, is collected with a separate trigger without any requirements on 
Hj or #x- The muon transverse momentum threshold is raised to Pt(f) > 35 GeV, while the 
Hj threshold is lowered from 500 GeV to 200 GeV and the jet multiplicity requirement is re- 
versed, to be fewer than three jets. Given these requirements on Hj and on the jet multiplicity, 
this control sample is dominated by SM processes. It is found that the estimated background 
agrees well with the number of events seen in the signal region Lp < 0.15. 

6.3 Results of the L P method 

The Lp distributions in three ranges of Sj P , are displayed in Fig. ^ for muons (top) and elec- 
trons (bottom). Tables [6] and [7] list the numbers of events observed and the number of events 
expected from all SM processes as presented above, in the signal region, for the muon and 
electron channels, respectively. The predictions, along with the numbers of events observed in 
each range of Sy P and Hj, are also displayed graphically in Fig. [7] for muons and in Fig. [sjfor 
electrons. The uncertainties quoted in Table [7]correspond to the statistical uncertainty of the fit, 
while the predictions displayed in Fig. [H] include the total statistical and systematic uncertainty. 

All estimates of the total contribution expected from SM processes in the various bins in (Sj P ,Hj) 
are consistent with the numbers of events observed in the data, with no visible excess from a 
potential SUSY signal. The result is interpreted as a limit in SUSY parameter space in the con- 
text of the CMSSM in Section gj 
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Figure 6: Data and fit results for the predictions for the Lp distribution, for events in the search 
sample, in different Sj P regions. Top plots for the muon channel; bottom plots for the electron 
channel. Left: 250 < Sj P < 350 GeV, center: 350 < Sj P < 450 GeV, and right: Sj P > 450 GeV. 



Table 6: Event yields in data and MC simulation for the muon sample. The results in the 
columns labelled "Total MC" are listed for reference only. The corresponding uncertainties 
statistical only. 



Sj P range 


Total MC 


Data 


Total MC 


SM estimate 


Data 


[GeV] 


Control region (Lp > 0.3) 


Signal region (Lp <0.15) 








500 < H T < 750 GeV 




[150-250) 


1465 ±11 


1297 


261 ± 3.2 


261 ±7 ±24 


258 


[250-350) 


452 ± 5.2 


383 


99.3 ±2.1 


84.1 ±4.2 ±7.3 


78 


[350^50) 


154 ±3.1 


128 


40.2 ±1.4 


33.3 ±3.0 ±2.6 


23 


>450 


59.2 ±1.8 


50 


18.6 ±1.0 


15.7 ±2.2 ±2.0 


16 






750 < Hj < 1000 GeV 




[150-250) 


280 ±4.1 


218 


52.4 ±1.6 


40.8 ±2.9 ±3.5 


46 


[250-350) 


91.9 ±2.1 


88 


22.3 ± 0.9 


21.3 ±2.3 ±2.2 


22 


[350^50) 


34.6 ±1.3 


25 


10.3 ±0.6 


7.5 ±1.5 ±1.0 


8 


>450 


26.7 ±1.4 


18 


8.8 ±0.6 


5.9 ±1.4 ±0.7 


7 






1000 GeV < Hj 






[150-250) 


92.3 ±2.5 


76 


20.5 ±1.0 


16.9 ±1.9 ±1.7 


15 


[250-350) 


32.9 ±1.3 


31 


8.7 ±0.8 


8.2 ±1.5 ±1.0 


8 


[350^50) 


10.9 ±0.7 


7 


4.6 ±0.4 


2.9 ± 1.1 ± 0.6 


1 


>450 


11.9 ±0.8 


12 


4.6 ±0.5 


4.6 ±1.4 ±0.7 


2 
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Figure 7: Comparison of the number of events observed in the data and the expectations from 

the background estimation methods for the muon channel, in the different Sj P bins. Left: 
500 < Hj < 750 GeV; Center: 750 < H T < 1000 GeV; Right: H T > 1000 GeV. The error bars 
indicate the statistical uncertainty of the data only, while the green band indicates the total 
statistical and systematic uncertainty on the background estimate. 



Table 7: Event yields in data and predictions of the numbers of EWK and QCD events for the 
electron sample in bins of Hj. The sum of predicted EWK events and predicted QCD events in 
the control region is constrained to be equal to the total number of data events. The background 
estimate used in comparing to the yields in the data is the result of the procedure described 
earlier and is listed in the row labeled "SM estimate". The uncertainties for the QCD and 
EWK background estimates are statistical only. The uncertainties shown for the SM estimate 
are first the statistical uncertainty from the control region fit and second all other systematic 
uncertainties. 



Sj P range [GeV] QCD EWK Data 
Control region (L P >0.3) 


QCD 


EWK SM estimate 
Signal region (L P <0.15) 


Data 






500 < Hj < 750 GeV 








[150-250) 


184 ±33 


1122 ±45 


1306 


9.1 ±1.6 


170 ±7 


179±7±18 


204 


[250-350) 


66 ±15 


334 ±22 


400 


2.1 ±0.5 


63.3 ±4.1 


65.3 ±4.3 ±5.9 


71 


[350^50) 


26.6 ±7.6 


93 ±11 


120 


0.3 ±0.1 


19.2 ±2.3 


19.4 ±2.4 ±2.9 


29 


> 450 


17.1 ±5.1 


33.9 ± 6.6 


51 


0.2 ±0.0 


9.0 ±1.8 


9.2 ±1.9 ±1.7 


11 






750 


< Hj < 


1000 GeV 








[150-250) 


39 ±15 


210 ± 20 


249 


1.9 ±0.7 


35.1 ±3.3 


37.0 ±3.5 ±4.8 


37 


[250-350) 


5.8 ±5.5 


59.2 ±9.1 


65 


0.2 ±0.2 


11.0 ±1.7 


11.2 ±2.0 ±1.8 


13 


[350^50) 


0.0 ±0.0 


26.0 ±5.1 


26 





6.3 ±1.2 


6.3 ±1.2 ±1.5 


5 


> 450 


8.7 ±3.4 


22.3 ±5.0 


31 


0.1 ±0.0 


6.7 ±1.5 


6.8 ±1.6 ±1.5 


5 








1000 GeV 


< Hj 








[150-250) 


14.9 ±7.7 


62 ±10 


77 


0.7 ±0.4 


11.7±1.9 


12.5 ±2.2 ±2.4 


9 


[250-350) 


10.4 ±4.3 


20.6 ±5.4 


31 


0.3 ±0.1 


4.5 ±1.2 


4.8 ±1.5 ±1.1 


8 


[350^50) 


0.5 ±1.7 


11.5 ±3.7 


12 


0.0 ±0. 


2.6 ±0.8 


2.6 ±1.2 ±0.9 


1 


> 450 


4.4 ± 2.5 


6.6 ±2.9 


11 


0.0 ±0.0 


2.5 ±1.1 


2.6 ±1.3 ±0.9 


1 



7 The Artificial Neural Network method 
7.1 Overview of the method 

The Artificial Neural Network (ANN) method uses a multi-variate analysis to combine sev- 
eral event characteristics, other than into a single variable that distinguishes signal from 
background. Signal events then preferentially populate a signal region in the two-dimensional 
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Figure 8: Comparison of the number of events observed in the data and the expectations from 

the background estimation methods for the electron channel, in the different Sj P bins. Left: 
500 < Hj < 750 GeV; Center: 750 < H T < 1000 GeV; Right: H T > 1000 GeV. The error bars 
indicate the statistical uncertainty of the data only, while the green band indicates the total 
statistical and systematic uncertainty on the background estimate. 



plane of the ANN output (zann) an d and the sidebands in this plane provide an estimate 
of the residual background. 

Four input variables drive the ANN. The first two are n^ ets , the number of jets with pj > 40 GeV, 
and Hj, the scalar sum of the pj of each jet with pj > 40 GeV. The SUSY signal typically has 
heavy particles decaying via complex cascades, and as such, is likely to produce more jets and 
larger Hj than SM backgrounds. The third variable is A<p(ji,j2), the angle between the two 
leading pj jets in the transverse plane, which makes use of the greater likelihood that the two 
highest pj jets are produced back-to-back in SM than in SUSY events. The final variable is Mj, 
the transverse mass of the lepton and £t system. In tt and W+jets events, the lepton and tj 
generally arise from the decay of a W boson, and as a result, Mj peaks near the W boson mass, 
with larger values arising only when there are additional neutrinos from T or semileptonic 
decays. By contrast, in SUSY events, Mj tends to be greater than the W mass because of f*j due 
to undetected LSPs. 

Figure [9] shows the distributions of these variables for simulated SM and SUSY events. The 
most powerful input variable is Mj; ftjets and Hj also have considerable discriminating power. 
The A(j)()i,)2) variable is weaker, but it still improves the sensitivity of the search. Lepton pj also 
discriminates between the SM and SUSY, but it is not included in the ANN because its strong 
correlation with £x in the SM would spoil the background estimate. Additional variables either 
do little to improve sensitivity or introduce a correlation between Zann and J£x- The input 
variables have similar distributions in the muon and electron channels, so we choose to train 
the ANN on the two channels combined, and use the same ANN for both. In general, the SM 
simulation describes the data adequately apart from a possible small structure near 130 GeV in 
the Mj distribution. Reweighting the simulation to match the Mj distribution in data does not 
affect the results of the analysis. 

The ANN infrastructure uses standard ROOT utilities |42l . During training, weights are de- 
termined that minimize the root-mean-square deviation of background events from zero and 
signal events from unity. For the SUSY parameter space under study, our sensitivity depends 
only mildly on the details of the signal sample that trains the ANN. Specifically, for LM points 
through 13 [30], the sensitivity is comparable (less than 30% variation) whether the ANN is 
trained on LM0, LM6 or LM9, even though these three training samples have rather different 
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Figure 9: The distributions of Hj ets , Hj, Acp, and Mj for data (solid circles), simulated SM 
(stacked shaded histograms), LM3 (open circles), and LM6 (open triangles) events after pre- 
selection. The small plot beneath each distribution shows the ratio of data to simulated SM 
yields. The muon and electron channels are combined. 



22 



7 The Artificial Neural Network method 



CMS 



4.98 ftr 



Vs = 7 TeV 




ANN 



Figure 10: The Zann distribution of the data (solid circles) and simulated SM (stacked shaded 
histograms), LM3 (open circles), and LM6 (open triangles) events, after preselection. The small 
plot beneath shows the ratio of data to simulated SM yields. 

characteristics. We select LMO for training because it gives the best overall performance. The 
SM simulation provides the background sample. 



Figure 10 compares the distributions of Zann f° r data and SM simulation for all events sur- 



viving the preselection. The two distributions are consistent within the uncertainties. The SM 
contribution is concentrated at small values of Zann/ while the LM3 and LM6 SUSY distribu- 
tions, which are also shown, extend to high values of Zann where the SM is suppressed. 

We define two signal regions in the two-dimensional $j and Zann plane. One region, referred 
to as the "low-£x" signal region, has Zann > 0.4 and 350 < £x < 500 GeV, while the other, 
the "high-£x" signal region, has the same Zann range, but JSj > 500 GeV. The high-^x signal 
region minimizes the probability that the expected background fluctuates up to a LM6 signal 
when signal contamination is taken into account. We observe 10 events in the low-£j signal 
region and 1 event in the high-£j signal region. 

7.2 Background estimation using the ANN sidebands 

The sidebands in the two dimensional plane of £j and Zann provide a strategy for estimating 



the background. The signal and sideband regions are shown in Fig. 11 and are denoted A, 
B, C, and D for the low-£j signal region and A, B', C, and D' for the high-^j signal region. 
The choice of boundaries for the sideband regions balances the competing needs of statistics 
and insensitivity to signal contamination against preserving similar event compositions in the 
signal and sideband regions. 



The predicted yield in region D is given by 



N, 



D,pred 



Nb x N c 
N A ' 



(6) 



where N, is the yield in region i, and the predicted yield in region D' is defined similarly. This 
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Figure 11: The yields of simulated SM (left) and LM6 (right) events in the £j versus Zann plane. 
The regions D and D' are the low- £7 and high-^x signal regions. The sideband regions are also 
indicated. 
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Figure 12: The £x distributions of events in the Zann signal region (solid circles) and sideband 
(green bars) for simulated SM (left) and data (right) events. The distributions are normalized 
in the £x sideband, 150 < j£j < 350 GeV (regions A and C for the two distributions respec- 
tively). The rightmost histogram bin includes overflow. The small lower plots show the ratio 
of normalized sideband to signal yields. 

procedure is equivalent to using the £x distribution of the Zann sideband regions (A, B, and B') 
as a template for the £j distribution of events with high Zann (Q D and D'), normalized using 
the yields in regions A and C. We test this estimation procedure using SM simulation: Fig. 12 
(left) shows that the jSj distributions for low and high Zann are similar. 

If a signal is present, it enters primarily in the signal regions D and D', but there are also signif- 
icant contributions relative to the SM in regions B and B', somewhat increasing the predicted 
backgrounds in D and D'. This effect is accounted for in the final results. 

Table [H] summarizes the event yields in the sideband subtraction regions for the various com- 
ponents of the SM background. The W+jets and tt dominate in all the regions, though their 
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Table 8: Event yields for the sideband (SB) and signal regions used in the ANN method. The 
uncertainties listed are statistical only. 

Sample type A B B 7 C D D 7 



tt single lepton 
tt dilepton (£) 
tt dilepton (t -+ £ 
tt single t 
Total tt 
W+jets 

Single top quark 
Z+jets 

QCD multijet 
Total SM 
Data 
LM3 

LM6 

relative proportion varies. The W+jets events are most important at low Zann since Mj, which 
largely drives zann/ tends to peak near the W-boson mass. Because the W bosons (and hence 
their daughters) can be highly boosted, these events extend to very high values of £x- As seen 
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zann signal 


zann signal 
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Low # T 
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Z.1V o 




u.z u.z 




17 4-07 
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u - u -o.o 
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o o+ - 2 
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o o +0 - 2 

u - u -o.o 
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5.4 ±1.2 


0.5 ±0.3 


171 ±7 


5.3 ±1.2 


0.4 ±0.3 


186 ±3 


20.4 ±1.1 


5.8 ±0.6 


40 ±2 


4.1 ±0.5 


1.6 ±0.3 


20 ±1 


1.5 ±0.3 


0.2 ±0.1 
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0.9 ±0.2 


0.1 ±0.1 


2.1 ±0.3 


07 +a12 
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0.8 ±0.1 
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3 +0 - 4 
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0.1 ±0.1 


oo +0 - 02 
u - uu -o.oo 
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u - uu -o.oo 
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224 ±7 


10.3 ±1.3 


2.1 ±0.4 


433 
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1 


164 ±3 


21 ±1 
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44.6 ± 0.5 


32.1 ± 0.4 


21.0 ±0.3 



in Fig. 10 tt events are more likely to have high values of Zann than are W+jets events; this 
is because of the presence of dilepton tt events, in which both W bosons (from the top quark 
pair) decay leptonically, but only one lepton is identified (dilepton (£)), giving large Mj. There 
is also a small contribution from events in which the lepton comes from the decay of a T pro- 
duced from a top quark decay, with the other top quark decaying either leptonically (dilepton 
(t — > £)) or hadronically (single t). The remaining small backgrounds come from single-top- 
quark, QCD multijet and Z+jets events. 

There are too few events in the simulated QCD multijet and Z+jets samples to populate the 
high regions (B, B', D and D'). For the results quoted in Table [8]for QCD multijet and Z+jets 
events, we employ an extrapolation technique based on loosening the Zann and require- 
ments. The extrapolated numbers for all the regions are consistent with those obtained from 
the simulated samples. The simulated yields in the sideband and signal regions indicate that 
QCD multijet and Z+jets events are negligible. 

The total SM simulation yields agree well with data in all regions, suggesting that the data 



share the main features described above. The Zann an d £t distributions are shown in Fig. 13 
7.3 Results of the ANN method 



Figure 12 (left) shows the results of applying the background estimation method to the SM 
simulation. We find that the method correctly predicts the background within a factor of k = 
D/ / D pred of °- 82 ± °- 12 ( stat -) in the low-£ T signal region and 0.69 ± 0.16 (stat.) in the high-£ T 
signal region. The modest deviation from unity results from a correlation between Zann and 
£t that arises because the W+jets background, which extends to large £t values, dominates in 
the Zann sideband (because it tends to have Mj near the W mass), whereas dileptonic tf events, 
with their somewhat softer £j spectrum, dominate in the Zann signal region. 



Figure 12 (right) shows the £t distributions of the data in the high and low Zann regions, after 



normalizing in the region 150 < £t < 350 GeV (A and C). Because the SM simulation appears 
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Figure 13: Distributions of j£j for (a) 0.2 < Zann < 0-4 and (b) z^nn > 0-4, and distributions 
of z ANN for (c) 150 < j£j < 350 GeV and (d) > 350 GeV. The samples shown are data 
(solid circles), simulated SM (stacked shaded histograms), LM3 (open circles), and LM6 (open 
triangles) events. The small plot beneath each distribution shows the ratio of data to simulated 
SM yields. 
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8 Systematic uncertainties 



to describe the data well, with, for example, consistent exponential decay constants describing 
the £x distributions in the ANN sidebands, we choose to scale the background prediction of the 
data by k. The uncertainty in the background from the relative cross sections of SM processes 
and other effects is quantified in Section[8] In the low-£x signal region, we expect 9.5 ± 2.2 (stat.) 
events, and in the high-^x signal region 0.7 ± 0.5 (stat.) events. The observed yields are 10 and 
1 events, respectively, consistent with the background prediction. 



8 Systematic uncertainties 



Systematic uncertainties affect both the background estimates and the signal efficiencies. The 
sources of systematic uncertainty in the background predictions vary among the three methods, 
both because the final event selections differ and because the background estimation methods 
themselves differ. The systematic uncertainties stem from lack of perfect knowledge of the 
detector response and from uncertainties in the properties of the SM backgrounds. Common 



uncertainties for all methods are described in Section 8.1 while details that are specific to each 



method are given in Sections 8.2 8.3 and 8.4 for the LS , Lp, and ANN methods, respectively. 
Tables [9 10 and 11 list the main uncertainties associated with each method. The systematic 
uncertainties affecting the signal efficiency and luminosity, which are largely common to all 
methods, are described in Section 1831 



8.1 Common uncertainties in the background predictions 

Table 9: Sources of systematic uncertainties for the LS method and their effects on the back- 
ground prediction in bins of The full list of systematic uncertainties is given for Hj > 
750 GeV, and the total uncertainties are shown for Hj > 500 GeV and Hj > 1000 GeV. Each 
uncertainty is expressed as a change in the ratio of the predicted to the true number of events 
(evaluated with simulation). Uncertainties associated with the dilepton and QCD backgrounds 
are discussed in the text. The total uncertainty is the individual uncertainties summed in 
quadrature. 
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[450-550) 


>550 




(%) 


(%) 


(%) 


(%) 




H T > 750 GeV 








Jet and energy scale 


11 


13 


14 


16 


Lepton efficiency 


1 


1 


1 


1 


Lepton px scale 


1 


2 


6 


2 


cr(tt) and cr(W) 


1 


1 


4 


4 


W polarization in tt 


1 


1 


1 


1 


W polarization in W+jets 


3 


4 


12 


11 


Z+jets background 


4 


4 


4 


4 


SM simulation statistics (K-factors) 4 


7 


12 


17 


Total systematic uncertainty 


13 


16 


24 


27 




H T > 500 GeV 








Total systematic uncertainty 


16 


18 


29 


30 




Hj > 1000 GeV 








Total systematic uncertainty 


15 


18 


28 


32 



8.1 Common uncertainties in the background predictions 



27 



Table 10: Sources of systematic uncertainties for the Lp method and their effects on the back- 
ground prediction in bins of Sj P for the muon and electron channels. The full list of systematic 
uncertainties are given for the range 500 < Hj < 750 GeV, and the total uncertainties are 
shown for the two ranges 750 < Hj < 1000 GeV and Hj > 1000 GeV. The total uncertainty is 
the individual uncertainties summed in quadrature. 



Sj P range [GeV]: 


[150-250) 


[250-350) 


[350-450) 


>450 




(%) 




(%) 




(%) 




(%) 




Channel 


¥ 


e 


¥ 


e 


¥ 


e 


¥ 


e 




500 < H T < 750 GeV 












Jet and £t energy scale 


6 


6 


4 


5 


5 


9 


9 


9 


Lepton efficiency 


5 


5 


5 


2 


3 


1 


1 


2 


Lepton pj scale 







1 




1 




2 




cr(tt) and cr(W) 


3 


1 


1 


1 


1 


2 


1 


1 


W polarization in tt 





1 


1 


1 


1 


1 


1 


2 


W polarization in W+jets 


2 


1 


2 


1 


2 


3 


3 


4 


i?x resolution 


2 


2 


1 


1 


1 


2 


4 


4 


tt(^) 


5 


5 


5 


5 


3 


3 


1 


1 


SM simulation statistics 


1 


1 


2 


2 


4 


5 


6 


7 


Total systematic uncertainty 


11 


10 


9 


8 


8 


12 


13 


13 




750 < Hj < 1000 GeV 










Total systematic uncertainty 


9 


12 


10 


11 


13 


13 


12 


13 




Hj > 1000 GeV 












Total systematic uncertainty 


10 


15 


13 


15 


20 


18 


16 


20 



The jet energy scale (JES) and its effect on $j in the event can affect the Hj and %j distributions 
and can also lead to differences between the lepton pj spectrum and f>j spectrum. To under- 
stand the effects of energy-scale variations, we vary the jet energy scale as a function of pj and 
rj by amounts determined in independent studies of jet energy scale uncertainties [38 |, and cor- 
responding to 2 GeV or less for jets with pj > 40 GeV, and then recompute Hj and £x- We also 
vary the energy scale of "unclustered" calorimeter deposits by 10% to take into account energy 
not clustered into jets (this effect is very small). 

The uncertainty in the lepton efficiency accounts for differences between data and simulation 
and uncertainties in the trigger efficiencies. The lepton efficiencies are studied using a sam- 
ple of lepton pairs with invariant mass close to the Z peak, in which one lepton satisfies tight 
selection criteria, and the second, reconstructed with relaxed criteria, serves as a probe of the 
tighter reconstruction and isolation requirements ("tag-and-probe" method [43 j). Discrepan- 
cies between the data and simulation for electrons are maximal at low pj (10% effect at around 
20 GeV), and we reweight events as a function of lepton pj to quantify the effect. The total 
lepton efficiency in data is described by simulation with an accuracy of 3%. Studies of the trig- 
ger that separately determine the efficiencies of the Hj 18ser , |? T tn 8S er , an d lepton requirements 
show that the lepton inefficiencies dominate, and amount to 2% to 3% for leptons that are re- 
constructed successfully offline. Muon pj scale uncertainties are obtained from the study of 
the q/pj (transverse curvature with sign given by the electric charge q) distribution of muons 
in Z events in data. By comparing the q/pj distribution of positive and negative muons it is 
possible to quantify the amount of bias in the measurement of q/ pj. 

The relative amount of tf and W+jets background affects each analysis method through correc- 
tions obtained from simulation. The contributions from tf and W+jet have not been specifically 
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8 Systematic uncertainties 



Table 11: Sources of systematic uncertainties for the ANN method and their effects on the back- 
ground prediction in bins of J&r- The total uncertainty is the individual uncertainties summed 
in quadrature. 



Ft rancrp TCtpVI' 

9-* I ACl-LL^C [VJC V J . 


r350-500) 


>500 




(%) 


(%) 


Jet and £x energy scale 


3 


4 


Lepton pi scale 


3 


5 


Lepton efficiency 


0.3 


0.4 


cr(tt) and cr(W) 


3 


2 


W polarization in W+jets 


1 


3 


W boson pj spectrum in W+jets 


10 


2 




1 


7 


Other backgrounds 


1 


1 


SM simulation statistics 


15 


23 


Total systematic uncertainty 


19 


26 



measured in the narrow region of phase space studied in this analysis and their relative con- 
tribution must be evaluated. The tf cross section is validated using an algorithm based on the 
reconstructed top-quark masses for both the hadronic and the leptonic top-quark decays. The 
uncertainty in the tt cross section is determined by comparing yields in data and simulation af- 
ter a selection based on top mass variables. The W+jets cross section is validated by comparing 
event yields between data and simulation in Z+jets events in a dedicated dilepton event selec- 
tion with similar kinematics. We assign an uncertainty to the W+jets cross section based on the 
agreement of the data and simulation in the Z+jets sample. Using the uncertainties obtained 
for the tt and W+jets cross sections, we probe different relative contributions of tt and W+jets 
events in our sample and the effect on our background predictions. 

Uncertainties in the polarization fraction for the W boson, either in tt or W+jets events, must 
be taken into account. For the W polarization in tt events, the theoretical uncertainties are very 
small (see Section|3]l and have negligible effect on the background predictions. The W polariza- 
tion in W+jets events, which is described in more detail in Section |3j is more complicated than 
in tt production. In this case, we consider the effect of conservative variations of the helicity 
fractions in bins of W-boson pj and rj with respect to the theoretical NLO calculations l|36| . 

For the dilepton tt background, tl(££), the uncertainties are evaluated somewhat differently for 
the different methods. In the Lp and ANN methods this background is evaluated together with 
the same control sample as for the main single-lepton background prediction. Uncertainties 
in the prediction can arise from finite detector acceptance, inefficient lepton identification, and 
cross section uncertainties. In the LS method the dilepton tf background is not predicted using 
the single-lepton background prediction and separate control samples must be used. Thus the 
uncertainties for the dilepton tt background are estimated separately and described in the next 
section. 

The small residual QCD multijet background is probed by inverting the requirement on I^°j mb 
or the electron selection criteria to obtain QCD dominated control samples. Contamination 
from tt and W+jet events in these control samples must be considered and the uncertainties on 
their cross sections are the dominant uncertainty for these methods. 

The Z+jets contribution to the signal regions is very small and uncertainties on this background 
prediction come from lepton efficiency and cross section uncertainties. In addition, for the LS 
method there is a small Z+jets contamination to the single-lepton control sample, which must 



8.2 Lepton Spectrum method background prediction uncertainty 
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be subtracted, and lepton efficiency and cross-section uncertainties are considered for this as 
well. 

8.2 Lepton Spectrum method background prediction uncertainty 

For the LS method, the systematic uncertainties for each of the different background predic- 
tions from control samples in data (1 I, dilepton, 1 t, QCD, and Z+jets) are included in Ta- 
bles [2j[3j and|4j To determine the systematic uncertainties for the largest source of background, 
l-£ events (arising from tf, W+jets, and single-top processes), we evaluate deviations for the 
^independent correction factor, which is determined from simulation and applied to the \-l 



background prediction (see Section 5.2 1. Table|9]gives a breakdown of the contributions of the 
systematic uncertainties for the l-£ prediction in bins of fij and for Hj > 750 GeV. The uncer- 
tainties in the l-l prediction for the Hj > 500 GeV and Hj > 1 TeV signal regions are similar 
to those listed in Table [9] The largest source of uncertainty arises from the potential difference 
in the muon pj and the £j scales, because the muon pj spectrum is used to predict the $j spec- 
trum. The statistical uncertainties in the correction factors (denoted as K-factors in Table [9} for 
the \-i method are slightly smaller than the combined systematic uncertainty of the correction 
factor. Table [9] does not include an uncertainty from jet resolution effects because this is taken 
into account by the smearing of the lepton pj spectra by QCD multijet £x templates (described 
in Section [5T2] |. For the purposes of setting limits, the total systematic uncertainty in the l-£ 
background prediction is treated as correlated across all bins in £x- 

Tables |2j|3j and [4] also list the non-single-lepton backgrounds, which account for about 25% of 
the total, with a relative uncertainty of 5-10% in the lowest- bin and about 30% in the highest- 
bin. For the dilepton prediction of lost and ignored leptons (described in Section |5.3| the 
main sources of systematic uncertainty arise from the lepton reconstruction and identification 
efficiencies and the top-quark pj spectrum. The uncertainties on the lepton efficiencies are 
described in Section |8.1[ and the uncertainty associated with the top-quark pj spectrum is 
determined from varying the fraction of events in the tail of this distribution in simulation 
in a manner consistent with the uncertainty in this tail as observed in data. This uncertainty is 
then propagated through the background determination procedure. 

8.3 Lepton Projection method background prediction uncertainty 

For the Lp method, the estimate of the total number of events expected from SM processes in the 
signal region, Ng^ d (Lp < 0.15), relies on the knowledge of the translation factor, Res/ as we U 
as the number of events observed in the control region, subtracted for the QCD background, 
Ndata(^p > 0.3). There are, therefore, two sources of uncertainty in this estimate: uncertainties 
in the number of events from EWK processes in the control region and uncertainty in Res- 
The relative change on the predicted background from each source of systematic uncertainty 

is listed in Table 10 for both muons and electrons. The largest uncertainty for high Sj P bins 
is the statistical uncertainty in the data in the control region. The second largest uncertainty 
comes from the JES uncertainty. The effect from the JES uncertainty is larger in the electron 
channel, since the JES affects also the shape of the Lp distribution used in the fit of the control 
region. The hadronic-recoil resolution uncertainty is conservatively evaluated by smearing the 
hadronic recoil by 7.5% in the parallel and 3.75% in the orthogonal direction. This decreases 
the resolution more than 10% for the high recoils (above 250 GeV) of the signal region and thus 
covers the difference between data and simulation. 
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9 Results and interpretation 



8.4 ANN method background prediction uncertainty 

For the ANN method, the systematic uncertainty in the background prediction is dominated 
by the statistics of the simulation, which probes for bias in the background estimation. Another 
important uncertainty comes from the pj spectrum of the W boson in W+jets events, since it 
affects the £j distribution of these events, which preferentially populate the zann sideband. 
To assess the impact, we reweight the pj spectrum of W boson events, using the differences 
in the pj spectra of Z bosons in data and simulation as a guide. This uncertainty is driven by 
the statistics of the Z+jets sample. The relative proportions of W+jets and tt events differ in the 
Zann signal and sideband regions so the background prediction depends on their relative cross 
sections. Those tt events with two leptons in the final state, only one of which is observed, have 
large #t and are the source of most SM events in the signal region. In addition to the tt cross 
section, this background depends on lepton acceptance and identification inefficiencies. Addi- 



tional sources of systematic uncertainty are the hadronic and leptonic energy scales. Table 11 
summarizes these uncertainties. 

8.5 Signal efficiency and other multiplicative uncertainties 

The systematic uncertainty in the signal yield arises from the uncertainty in the signal efficiency. 
In general, this uncertainty is correlated across Ifcj or Sj P bins. The JES component of the signal 
efficiency uncertainty is computed separately for each model point in CMSSM and simplified 
model parameter space and is correlated with the JES uncertainty in the single-lepton back- 
ground prediction. The systematic uncertainties in the signal efficiency associated with lepton 
reconstruction and the trigger amount to 3%. The uncertainty in the integrated luminosity is 
2.2% [44] . The systematic uncertainty in the signal efficiency, not including the JES component, 
is 6% for each of the analyses. 



9 Results and interpretation 

The LS, Lp, and ANN methods each yield SM background predictions that are consistent with 
the number of events observed in data. We therefore proceed to set exclusion limits on SUSY 
model parameters. All limits are computed using the modified-frequentist CLs method l|45| 
with a one-sided profile likelihood test statistic. To interpret the absence of an observed signal, 
three complementary approaches are used. 

9.1 Constraints on CMSSM parameter space 

First, we scan over models in the CMSSM and determine whether the number of events pre- 
dicted at each model point in parameter space can be excluded by the measurements. This 
procedure relies on the fact that the CMSSM parameter space can be described with just five 
parameters, and we fix three of them to commonly used values (Aq = GeV, \i > 0, tan f> = 10). 
Each model point has a complete SUSY particle spectrum and a well defined cross section, 
which typically involves several production subprocesses. The CMSSM simulated samples 
are initially generated using leading-order cross sections. At each point in CMSSM parame- 
ter space, the predicted yields for each production subprocess (e.g., gg — > gg) are corrected 
using the NLO cross sections discussed in Ref. H46ll . Using the observed yield in data and the 
predicted background, we determine whether the CMSSM yield for the particular model point 
can be excluded at the 95% confidence level (CL). This procedure is complicated by the fact 
that the control regions in data could potentially be contaminated by signal events. This effect 
is taken into account for each model by removing the expected contribution to the predicted 
background arising from signal contamination of the control regions. 



9.2 Constraints on simplified model parameter space 
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Figure 14: LS method: exclusion region in CMSSM parameter space for the Hj > 750 GeV 
selection. 



Figures 14 15 and 16 show the CMSSM exclusion region ||47| for the three background esti- 
mation methods, evaluated in the mi/2 vs. mo plane, with the values of the remaining CMSSM 
parameters fixed at tan ft = 10, Aq = GeV, and }i > 0. Figure [17] displays all of the results 
together. The excluded regions are below the plotted curves, corresponding to SUSY parti- 
cle masses below certain values. For reference, the plots display curves of constant gluino 
and squark masses. The lines of constant gluino mass are approximately horizontal with 
m (g) 2.5 trii/2- Lines of constant squark mass are strongly curved in the mi/ 2 vs. niQ plane. At 
low wiq, the analyses exclude gluinos with masses up to about 1.3 TeV, but the sensitivity falls 
with increasing niQ. To determine the one standard deviation (a) theoretical uncertainty on the 
observed limit, the signal yields are recomputed after changing each of the process-dependent 
SUSY production cross sections at each model point by ±1(7 of their uncertainty arising from 
the parton distribution functions and renormalization and factorization scales |46|. 

9.2 Constraints on simplified model parameter space 

The second approach to interpretation is based on the use of simplified models |T4l[T5| , which 
provide a more generic description of new physics signatures. Such models do not include a 
full SUSY particle spectrum, but instead include only the states needed to describe a particular 
set of decay chains of interest. Rather than excluding a model, the procedure is to calculate cross 
section upper limits on a given topological signature. (Such cross section limits can, however, 
be converted into limits on particle masses within the assumptions of the particular model.) 
Because simplified models do not describe a full SUSY spectrum, the number of free parameters 
is small. Furthermore, the parameters are simply the masses of the SUSY particles, in contrast 
to the grand-unified-theory-scale parameters used in the CMSSM. An advantage of simplified 
models is that, as a consequence, certain relationships between particle masses that arise with 
the CMSSM no longer hold, and the spectra can be much more generic. 
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Figure 15: Lp method: exclusion region in CMSSM parameter space for all Hj bins combined. 




We consider the "Topology 3 weakino" (T3w) simplified model, which involves the produc- 



tion of two gluinos and their decay via the mechanism shown in Fig. 18 One gluino is forced 
to decay into two quark jets plus the LSP (x°) via the three-body decay g — > qqx°, while the 
other gluino decays via g — > qq'^, followed by x — > W ± x°- The W ± boson can then de- 
cay leptonically. The T3w model is specified by masses of the gluino, the LSP (x°), and an 
intermediate chargino (x ). We calculate cross section limits as a function of M(g), assuming 
a fixed value for the LSP mass M(x°) = 50 GeV and setting the chargino mass according to 
M(x ± ) = 0.5(M(x°) + M(g)). The nominal production cross section for the gluino pair pro- 



duction mechanism is given in Ref. |46| . Figure 19 shows the cross sections excluded by each 
method for this model. The limits fluctuate significantly at low M(g) because of the low signal 
efficiency in this region. 

9.3 Alternate model exclusions 

The data can be interpreted using a third approach, which is applicable to models that do not 
fall within the scope of either the CMSSM or the simplified model discussed in this section. 
A model builder can investigate the sensitivity of any one of the three methods presented in 
this paper to a given signal hypothesis by applying the event selection requirements listed in 
Table [l] together with the final requirements that define the signal regions. We provide a simple 
efficiency model for the most important observables used in the event selections. The efficiency 
model can then be applied to a basic (pythia) simulation of the signal process. 

The efficiency model is based on parametrizations of the efficiencies for the event selection re- 
quirements with respect to the main reconstruction objects and quantities, such as Hj, and 
lepton pj. The efficiency of the analysis for a given model can be estimated by applying these 
individual reconstruction efficiencies, which are given as a function of the most important pa- 
rameter (such as lepton pj), to the corresponding kinematic distributions in the model. This 
procedure would then yield an estimate for the number of signal events from the model. Fi- 
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Figure 16: ANN method: exclusion region in CMSSM parameter space. 
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Figure 17: Exclusion region for the LS, Lp, and ANN methods in CMSSM parameter space. 
Results from the low- and high-^x signal regions are combined. 
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Figure 19: Excluded cross sections for the LS, Lp, and ANN methods for the T3w simplified 
model. 
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nally the sensitivity of the analysis to the model can be obtained by comparing the yield of sig- 
nal events obtained in this manner with the background yields given in this paper. Kinematic 
correlations (which can be model dependent) are not taken into account, but this approach 
nonetheless provides a first approximation to the sensitivity. 

The efficiencies for each analysis object are described using "turn-on" curves, which are simply 
error functions, 



^plateau 2 



erf( * + ! 



(7) 



where x represents the variable most relevant for the reconstruction of the particular object. The 
error function is parametrized in terms of the plateau efficiency, e p i a teau/ the turn-on threshold, 
^thresh; and the characteristic width of the turn-on region, a. These parameters are obtained by 
fitting simulated event samples as a function of the true (generated) value. 

The selection efficiency associated with the lepton reconstruction, identification, and isolation 
requirements is estimated as a function of lepton pj by considering muons and electrons (in- 
cluding those from t decay) generated in the PYTHlA-simulated hard-scattering process. The 
lepton isolation requirement has a large effect on the efficiency, which consequently depends 
on the number of jets in the event. To reduce the model dependence arising from this effect, 
two categories of leptons are considered. First, we assign zero efficiency to leptons that are 
within AR < 0.4 of a quark or gluon with pj > 40 GeV in the hard-scattering process. The ef- 
ficiency for the remaining leptons is described by a turn-on curve whose parameters are listed 
in Table (12) The efficiencies are specified for both the lepton selection and for the lepton veto. 

Table 12: Efficiency-model parameters for lepton efficiencies as a function of x = pj. The 
leptons are required to lie within the fiducial region and must satisfy the pj thresholds specified 
in Table [I] 



Lepton 


^plateau 


^thresh [GeV] 


cr [GeV] 


Muon (signal) 


0.86 


2.7 


65 


Muon (veto) 


0.90 


-17 


75 


Electron (signal) 


0.74 


20 


61 


Electron (veto) 


0.83 


2.3 


54 



The number of jets and the resulting Hj value for each event are computed using information 
available at the generator level. The same clustering algorithm used to reconstruct jets in the 
data is applied to the generator-level particles. The resulting generator-level jets are required 
to satisfy AR > 0.3 with respect to the leptons described above. The £t variable is estimated at 
the generator level from the transverse momenta of neutrinos and any new weakly interacting 
particles, such as the The parametrizations of the efficiency turn-on curves for the Hj and 



£t requirements are listed in Tables 13 and 14 respectively. For the requirements used with the 
LS method, the information given in these tables generally reproduces the efficiency from full 
simulation to within about 15%. 

In the Lp method, the variables Lp and Sj P are functions of lepton pj and J£j- The modeling 
of lepton pj is described above. To emulate #t, one needs to apply both a scale shift and 
smearing to the generated $j value. The J£j scale factor is £x( r eco)/^x(gen) = 0.94. The 
value of cr(£ T (reco)/£ T (gen)) is about 0.2 at £ T = 100 GeV. It falls linearly to about 0.06 at 
It = 400 GeV, and it remains at 0.06 for j£ T > 400 GeV. 

In the ANN method, the preselection requirements on Hj and £x are 400 and 100 GeV, respec- 
tively. The signal regions are specified by 350 < j£j < 500 GeV and £j > 500 GeV together with 
Zann > 0-4, where Zann is a function [48] of rij ets , Hj, A(p()i,)2), and Mj. The efficiency turn-on 
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curve for zann > 0.4 is approximated by the parameter values e p i a teau = 0.98, x t hresh = 0.41, 
and a = 0.1. 

With these additional procedures, the emulation of the efficiencies for the Lp and ANN meth- 
ods is found to be accurate to within ~15%, as for the LS method. 



Table 13: Efficiency-model parameters for x = H T . 



Threshold 


^plateau 


^thresh [GeV] 


a [GeV] 


Hj > 400 GeV 


1.00 


396 


65 


Hj > 500 GeV 


1.00 


502 


66 


Hj > 750 GeV 


1.00 


760 


68 


Hj > 1000 GeV 


1.00 


1013 


80 


Table 14: Efficiency-model parameters for x = %j. 


Threshold 


^plateau 


Xthresh [GeV] 


a [GeV] 


£ T > 100 GeV 


1.00 


103 


41 


£ T > 250 GeV 


0.99 


266 


41 


£ T > 350 GeV 


0.98 


375 


45 


£ T > 450 GeV 


0.97 


485 


48 


£ T > 500 GeV 


0.94 


537 


44 


£ T > 550 GeV 


0.96 


597 


59 



10 Summary 

Using a sample of proton-proton collisions at yfs = 7 TeV corresponding to an integrated lumi- 
nosity of 4.98 fb _1 , we have performed a search for an excess of events with a single, isolated 
high-px lepton, at least three jets, and large missing transverse momentum. To provide a robust 
and redundant determination of the SM backgrounds, three methods are used, each of which 
relies primarily on control samples in the data. 

The Lepton Spectrum (LS) method exploits the relationship between two key observables, the 
lepton px distribution and the £x distribution. In the dominant SM background processes, 
which have a single, isolated lepton, this connection arises from the fact that the lepton and 
neutrino are produced together in the two-body decay of the W boson, regardless of whether 
the W is produced in tt or W+jets events. In many SUSY models, however, the £j is associated 
with the production of two neutralinos, which decouples £j from the lepton px spectrum. 
Smaller backgrounds arising from tt dilepton events, from x —> i decays in tf or W+jets events, 
and from QCD multijet processes are also estimated using control samples in the data. In the 
sample investigated with this method, at least four jets are required, which helps to suppress 
the background from W+jets events. Nine signal regions are considered, specified by three 
thresholds on Hj and three bins of %j. The observed yields in each region are consistent with 
the background estimates based on control samples in the data. 

The Lepton Projection (Lp) method exploits information on the W-boson polarization in tt and 
W+jets events. The dimensionless Lp variable itself is sensitive to the helicity angle of the lepton 
from W decay, but it also provides discrimination between signal and background through 

the ratio of the lepton px and the £x values, which is small in SUSY-like events. The Sj P 
variable maps out a diagonal line in the plane of lepton p T vs. £t and reflects the W transverse 

momentum for the boosted W boson. The Lp distributions are studied in bins of Sj P , and 
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Hj, and at least three jets are required. In each signal region, the data are in agreement with 
expectations from the SM. 

The artificial neural network (ANN) method provides a means to obtain the £t distribution of 
background events in data by constructing a neural network variable ZanN/ which has a very 
small correlation with j£x- This variable also provides strong discrimination between signal and 
background events, so that the background regions do not suffer from large signal contamina- 
tion in the models considered. A key element of the Zann variable is the transverse mass of the 
lepton-£j system, but additional variables, such as the number of observed jets, play a role as 
well. In the ANN analysis, no excess of events is observed in the signal regions with respect to 
the SM background prediction. 

The results from each method are interpreted in the context of both the CMSSM and a so- 
called simplified model, T3w, which has a minimal SUSY particle spectrum. The CMSSM limits 
exclude gluino masses up to approximately 1.3 TeV in the part of the parameter space in which 
mo < 800 GeV, but the bound gradually weakens for larger values of mo- For the T3w simplified 
model, we obtain cross section upper limits as a function of gluino mass. Finally, we provide 
an approximate model of our signal efficiency that can be used in conjunction with a simple 
PYTHIA simulation to determine whether other models can be probed by these data. 
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